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Our  concern  is  with  the  flow  around  an  infinite  cylinder,  which  at  a  certain  instant  is 
impulsively  started  to  spin.  The  growth  of  vortices  in  the  resulting  boundary  layer  occurring 
outside  the  cylinder  is  investigated.  This  layer  is  essentially  a  Rayleigh  layer  which  grows 
with  time,  so  the  mechanism  involved  is  similar  to  that  studied  in  Hall  (1983).  Vortices 
with  wavenumber  comparable  to  the  layer  thickness  are  shown  to  be  described  by  partial 
differential  equations.  It  is  found  necessary  to  solve  the  unsteady  partial  differential  equations 
that  govern  the  system  numerically.  We  assume  that  the  Rayleigh  layer  is  thin  so  particles 
are  confined  to  move  in  a  path  with  radius  of  curvature  the  same  as  the  cylinder.  The  Gortler 
number  is  a  function  of  time,  so  we  consider  the  time  scale  which  produces  an  order  one 
Gortler  number.  We  consider  the  right  hand  branch  calculation  by  letting  the  time  tend  to 
infinity,  also  inviscid  Gortler  modes  are  considered. 
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§1  Introduction 


The  problem  we  are  considering,  is  that  of  an  infinite  cylinder  held  fixed,  until  at  a 
certain  instant  when  a  torsional  impulsive  velocity  is  imparted  to  the  cylinder,  without 
loss  of  generality  this  time  may  be  considered  to  be  t  =  0.  The  cylinder  is  immersed 
in  a  fluid  of  kinematic  viscosity  v  and  density  p,  and  the  cylinder  has  radius  a.  The 
torsional  velocity  imparted  to  the  cylinder  is  Cl.  It  is  known  that  a  layer  of  thickness 
a(vt)*  occurs  at  the  cylinder,  this  is  a  Rayleigh  layer  which  was  shown  to  occur  on 
a  infinite  flat  plate  by  Stokes  (1855)  and  later  by  Rayleigh  (1911).  We  consider  small 
time  scales  so  that  the  resulting  layer  is  thin  and  thus  the  particles  contained  in  this 
layer  can  be  considered  to  move  in  paths  with  the  same  radius  of  curvature  as  the 
cylinder.  We  may  define  the  Gortler  number  as 


in  the  same  way  as  in  Seminara  &  Hall  (1976)  and  Hall  (1984).  These  papers  considered 
centrifugal  instabilities  in  Stokes  layers,  so  (i 4) 5  is  replaced  by  (u/u)^ ,  where  u>  is  the 
frequency  of  the  oscillations.  The  Reynolds  number  is  defined  here  by 


(1.2) 


The  ensuing  problem  is  similar  initially  to  Hall  (1983),  who  considered  a  spatially 
growing  boundary  layer  occurring  on  a  concave  wall.  We  wish  to  investigate  vortices 
similar  to  those  found  to  occur  by  Taylor  (1923)  and  Gortler  (1940).  In  Hall  (1982 
a,b)  a  self-consistent  asymptotic  description  of  the  growth  of  short  wavelength  Gortler 
vortices  was  given.  In  particular  in  Hall  (1982  a),  it  was  shown  that  the  equations  gov¬ 
erning  the  centrifugal  instabilities  are  partial  differential  and  not  ordinary  differential 
equations  as  discussed  by  Gortler  (1940)  and  subsequent  authors.  In  Hall  (1982  b)  the 
nonlinear  development  of  the  vortices  was  considered.  In  this  paper  we  concentrate  on 
disturbances  with  wavelengths  initially  of  order  one,  but  will  consider  the  right  hand 
branch  in  terms  of  local  wavelength  and  local  Gortler  number. 
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We  may  now  apply  the  usual  non-dimensionalisation  and  non-dimensionahse  time 
with  respect  to  a2/ u,  thus  the  Gortler  number  in  (1.1)  may  be  written  in  terms  of  the 
Reynolds  number  Rei  and  the  non-dimensional  time  T  as 

G  =  R2eT%,  (1.3) 

To  invoke  the  mechanism  as  described  in  Seminara  &  Hall  (1976),  it  is  necessary 

_  4 

to  regard  G  as  order  one,  and  this  can  be  seen  to  occur  if  T  ~  i?e  3 ,  and  since  we 
consider  Re  >  1,  this  is  a  small  time  scale.  We  may  define  a  similarity  variable  77,  and 

thus  it  is  possible  to  solve  the  basic  flow  equations  in  the  resulting  layer,  note  that  we 

—  2 

consider  a  radial  scaling  of  Re  5  so  that 

f  =  =  0(1).  (1.4) 

where  R  is  a  dimensional  radial  variable.  The  essential  point  here  is  that  the  equations 
governing  the  motion  in  this  layer  involve  temporal  derivatives  as  well  as  radial  ones, 
this  is  equivalent  to  ensuring  ^-derivatives  are  retained  in  spatial  problems,  a  fact 
that  seems  to  have  been  first  noted  by  Smith  (1955).  Smith  (1955)  retained  some  of 
the  terms  associated  with  the  non-parallelism  of  the  Gortler  problem.  For  Tollmien- 
Schlichting  waves  Bouthier  (1973)  and  Gaster  (1974)  incorporated  extra  terms  to  ap¬ 
preciate  the  effect  of  boundary  layer  growth,  which  greatly  improved  the  comparison 
between  experimental  and  theoretical  results.  Although  it  was  noted  by  Gaster  (1974) 
that  the  Bouthier-  Gaster  work  just  amounts  to  a  successive  approximation  technique, 
Smith  (1979)  attempts  to  remedy  this  and  does  so  effectively.  Here  we  consider  an 
unsteady  instability,  this  is  included  in  Kerczek  &  Davis  (1974)  and  subsequently  in 
Kerczek  &  Davis  (1976)  and  Kerczek  (1982),  where  in  these  papers  they  consider  pe¬ 
riodic  flow.  The  description  of  the  flow  in  a  boundary  layer  growing  with  time  on  a 
semi-infinite  flat  plate  was  attempted  in  Otto  (1992  b),  employing  a  triple  deck  struc¬ 
ture  as  in  Smith  (1979  a),  so  it  is  known  theoretically  that  vortices  occur  in  a  time 
dependent  growing  boundary  layer. 

Here  we  concentrate  on  disturbances,  of  order  one  wavelengths  as  mentioned  pre¬ 
viously  as  these  are  most  likely  to  be  of  practical  importance.  The  development  of 


the  vortices  is  caused  by  the  inclusion  of  a  small  hump.  This  hump  is  considered  to 
be  the  same  size  as  the  boundary  layer,  this  technique  was  employed  in- Denier,  Hall, 
&  Seddougui  (1991)  in  the  spatially  varying  Gortler  problem.  Various  forms  of  this 
hump  are  chosen  and  give  rise  to  different  regions  of  vortex  activity,  subsequently  re¬ 
ferred  to  as  vortex  wedges.  In  this  paper  we  do  not  have  the  freedom  of  Denier,  Hall 
&  Seddougui  (1991)  to  introduce  the  disturbance  at  a  certain  spatial  position  as  we 
have  said  that  the  impulsive  motion  starts  at  time  equals  zero,  and  the  hump  must  be 
present  for  all  time. 

The  identification  of  the  time  station  T,  when  a  particular  vortex  wavenumber  is 
neutral  incurs  the  same  problem  as  in  non  parallel  spatial  calculations,  refer  to  Hall 
(1983),  Gaster  (1974)  and  Smith  (1979  a).  In  this  paper  the  flow  quantity  chosen  to 
be  representative  of  the  vortex  size  is 


CO 


0 

and  thus  the  neutral  position  is  defined  by  the  position  where  a,  defined  by 


(1.5) 


a  =  E 


(1.6) 


becomes  zero. 

Thus  at  this  time  station  we  may  calculate  the  local  wavenumber  a,f  and  the  local 
Gortler  number  Gf .  The  values  of  the  equivalent  pair  in  Denier,  Hall,  &  Seddougui 
(1991)  are  found  *o  depend  on  the  initial  form  of  the  disturbance.  The  initial  form  of 
the  disturbance  is  a  function  of  the  hump  at  the  wall. 

We  now  go  to  investigate  inviscid  Gortler  modes,  where  we  have  a.f  ~  1  and 
Gf  >  1  described  in  Denier,  Hall,  &  Seddougui  (1991).  We  also  consider  the  righthand 
branch  problem  such  that  for  (if  >  1,  we  have  Gf  ~  ci}F,  as  investigated  in  Hall  (1982 
a).  In  this  problem  we  have  constant  curvature,  thus  the  right  hand  branch  cannot  be 
encountered  by  letting  k  ~  X  ?  as  in  Hall  (1982  a).  In  order  to  investigate  the  right 
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hand  branch  it  is  found  necessary  to  consider  very  small  initial  wavelengths  k  at  a 
time  T  such  that, 

with  T»  1.  (1.7) 

We  also  consider  matching  between  the  inviscid  solutions  and  right  hand  branch  solu¬ 
tions  as  in  Denier,  Hall  &  Seddougui  (1991). 

The  approach  adopted  in  the  remainder  of  this  paper  is  as  follows;  in  section  2  we 
derive  the  linear  0(1)  wavenumber  regime  problem,  in  section  3  the  inviscid  Gortler 
problem  is  tackled,  section  4  contains  details  of  the  right  hand  branch  calculations, 
section  5  includes  description  of  the  numerical  techniques  involved  in  the  solution 
of  the  equations  derived  in  sections  3  and  4.  Finally  section  6  contains  results  and 
discussions  including  comments  on  matching  between  the  solutions  of  section  4  and 
section  5,  and  a  discussion  of  possible  further  developments  to  take  account  of  fully 
nonlinear  right  hand  branch  calculations.  A  discussion  of  extension  to  the  problem  to 
further  geometries  is  also  included. 
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§2  Formulation  of  linear  order  one  problem 

We  non-dimensionalise  the  velocity  field  with  respect  to  a£l,  and  the  pressure 
terms  with  respect  to  pa2Q?.  The  induced  motion  is  taken  to  be  independent  of  #,  the 
coordinate  system  employed  is  that  of  cylindrical  polars  with  components  (r,  9,  z ),  and 
corresponding  velocity  components  (u,  v,  w),  the  Navier-Stokes  equations  then  become, 

l  dU  rTdU  „rdU  V2  dP  1  ((d2  d2  \  Id  l\Tr  , 

Re  dt  8r  +W  dz  r  dr  +  Re  ((dr2  +  dz2)  +  r  dr  r2 )  U’ 

1  dV  dV  „rdV  VU  1  ( ( d2  d2\  Id 
TeW  +  U^  +  W^+~  =  Te{{d^  +  d^)+7^~^)^  (2,16) 


1  dW  TTdW  „rdW  dP  1 

+  U  +w  =-  + 

Re  dt  dr  dz  dz  Re 


dr 2  +  dz2 


+  1r£)W’ 


dU  U  dW  „ 

aT  +  7  +aT  _0- 


(2.1c) 


(2.1  d) 


We  now  introduce  an  axi-symmetric  hump  at  the  cylinder,  defined  by  r  =  1  +  6f  (z) 
.  Following  the  arguments  in  the  introduction,  it  is  found  necessary  to  look  at  a  time 
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scaling  T  ~  Rc  3  so  that  the  Gortler  number  is  order  one,  and  in  this  section  it  will  be 
taken  to  be  one  with  no  loss  of  generality.  Thus  we  define  boundary  layer  variables,  rj 
and  Z  given  by 

r  =  1  +  Rc3V,  z  =  Rc3Z ,  (2.2c) 

and  a  time  boundary  layer  variable  r,  defined  by 


T  =  Rc  3V, 


(2.2  b) 


To  determine  the  basic  flow  we  assume  it  to  be  independent  of  Z  and  a  function  of  £, 
a  similarity  variable  given  by 

<  =  A-  (2-3) 

2t  2 


The  basic  flow  is  found  to  be  governed  by 


d2v  dv 

d?  +2^“0, 


(2.4a) 
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with  boundary  conditions 


v  =  0  as  £  —*■  oo, 


v  =  1  at  £  =  0, 

The  solution  of  (2.4)  is  the  error  function,  given  by 


(2.46) 

(2.4c) 


OO 

5K)=^/e'?4 


We  now  introduce  a  disturbance  so  that  the  flow  field  is  now  given  by 

u  =  (0,u,0,p)  +  6(Re  au,v,Re  3w,Re  3p)  H - ,  (2.6) 

where  the  quantities  with  tildes  are  functions  of  rj,  Z  and  r,  and,  as  previously  stated 
independent  of  0.  We  take  6  to  be  a  small  vanishing  parameter,  so  that  we  are  in  a 
linear  regime,  with  this  assumption  the  governing  equations  are  found  to  be, 


dz  n__  op  ,  (  d2  ,  d2  \  _ 

dr  2VV  dp  +  [dp*  +  dZ 2  )  U ’ 

dv  dv_  ( d2  d2  \ 

dr  +  dt]U~  \dtf  +  dZ2)V’ 

dw  dp  f  d2  d2  \  _ 

Ik  =  ~dz  +  +  dz* )  Wy 

du  dw  n 
d^  +  dZ~0, 


(2.7  a) 


(2-7  b) 


(2.7c) 


(2.7  d) 


with  boundary  conditions  given  by 


u  =  v  =  u>  =  0  as  T)  —*  oo, 


(2.7  e) 


so  that  the  disturbance  is  confined  to  the  Rayleigh  layer, 


u  =  v  =  —  =  0  at  t]  =  f(Z). 
dr] 


(2-7 1) 


Here  (2.7  f)  represents  the  no  slip  boundary  conditions  at  the  hump  on  the  cylinder, 
which  is  chosen  to  be  scaled  with  the  boundary  layer  thickness  multiplied  by  the  dis¬ 
turbance  size  8.  Note  that  since  we  looking  close  to  the  cylinder  the  terms  proportional 


to  secondary  radial  corrections  are  neglected,  such  as  1/rJ^.  Let  us  now  apply  the 
Fourier  transform  to  the  disturbance  quantities,  thus  the  transform  U  of  u  is  defined 

by 

CO 

U  =  -=  J  ueikZdZ.  (2.8) 

—  CO 

Thus  the  system  (2.7)  now  becomes 


dU  ~  dP  (  d2  ,2\- 

dr  2vV  ~  dtj+  ( dr f  k  )  u' 


dU 

—  +  ikW  =  0, 
dt] 


(2.9a) 

(2.96) 
(2.9c) 
(2.9  d) 


with  boundary  conditions  of 


II 

'>- 

II 

'to 

W  =  0  as  rj  —y  oo, 

(2.9e) 

o" 

II 

II 

'to 

~  dv 

V  =  —F—  at  t]  =  0. 

drj 

(2.9/) 

where  F  is  the  Fourier  transform  of  the  hump  f(z).  Note  that  is  the  absence  of  the 
hump  there  would  be  no  vortex  mechanism.  Now  by  elementary  manipulation,  we 
may  eliminate  W  and  P,  to  obtain  the  coupled  partial  differential  system, 

(  d2  ,<>  d  \  (  d2  _ _ ^ 


(£__ 

w 


d 2 
dr? 

-lAu  =  VrvV, 

(2.10a) 

d\ 
dr  ) 

(2.106) 

this  system  is  the  Gortler  equations  for  a  vortex  with  temporal  growth,  with  Gortler 
number  equal  to  one.  It  is  now  noted  that  there  will  be  a  discontinuity  in  the  dis¬ 
turbance  flow  field  due  to  the  impulsive  start,  so  it  is  necessary  to  solve  the  problem 
asymptotically  for  small  r.  Since  the  solution  of  (2.10)  necessitates  a  numerical  ap¬ 
proach,  this  asymptotic  solution  may  be  used  as  an  initial  condition.  It  can  be  seen 
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J  w  3 

from  (2.10  b)  that  F(0,t)  ~  r~  2,  17(0,  r)  ~  r  j  for  r  <C  1  following  this  scaling  we 
seek  a  solution  of  the  form 


V  =  r  y0  (£)  +  ...  ,  (2.11a) 

U  =  T->U0(O  +  -"  ,  (2-116) 

where  £  is  the  similarity  variable  defined  by  (2.3),  it  may  be  shown  that  Vo(£)  satisfies 


d2V0 

d? 


dV0 

+  2^+2Vo=0, 


(2.12  a) 


with  boundary  conditions 


Vb(0)  =  2 


and  Vo(oo)  =  0, 


(2.126,  c) 


where  F(k)  is  the  Fourier  transform  of  f(z),  the  form  of  the  hump.  The  solution  of 
(2.12)  is  given  by 

TPflA 

(2.13) 


Note  that  we  take  the  radial  component  Vo  to  be  zero  for  small  t,  as  it  is  0(t2)  smaller 
than  Vo-  The  equations  may  now  be  solved  numerically,  the  solution  is  discussed  in 
section  5.  Note  that  since  this  system  is  linear  F(k )  may  be  taken  to  be  unity,  this 
corresponds  to  a  6-function  shaped  hump,  so  is  relevant  to  a  localised  imperfection  on 
the  surface  of  the  cylinder. 


§3  The  inviscid  limit  of  the  Gortler  problem 

In  Denier,  Hall  &  Seddougui  (1991)  the  inviscid  Gortler  problem  was  discussed, 
it  corresponds  to  a  large  Gortler  number  and  order  one  wavenumbers.  The  temporal 
inviscid  Gortler  problem  is  virtually  identical  to  the  inviscid  problem  set  as  an  exercise 
by  Drazin  &  Reid  (1979)  and  solved  by  Floryan  (1986).  Thus  we  have  the  governing 
equations  (2.10)  with  a  Gortler  number  other  than  one,  taken  to  be  G 


d2 

On2 


U  =  2  Gk2vV, 


(3.1  a) 


8 


a;;2  3i)  ’ 


(3.16) 


We  now  consider  a  solution  of  (3.1)  in  the  inviscid  regime,  that  is  G  1  with  k  =  0(1). 
It  may  be  shown  that  the  growth  rate  is  proportional  to  CH  following  Denier,  Hall  & 
Seddougui  (1991).  Thus  the  disturbance  in  this  regime  takes  the  form 


(U,V)  =  E(U,G-V), 


where 


E  =  exp  (  G 


*r'+ 


(3.2  a) 


(3-2  b) 


This  form  of  the  disturbance  may  be  substituted  into  (3.1),  and  the  two  equations  can 
be  combined  to  produce, 


**(J 


In  Denier,  Hall  &  Seddougui  (1991)  the  equivalent  equation  has  a  closed  form  of 
solution.  Unfortunately  (3.3)  does  not,  so  thus  must  be  solved  numerically.  The 
equation  (3.3)  was  solved  using  a  fourth  order  Runge-Kutta  scheme.  The  growth  rate 
is  plotted  against  k  in  figure  1.  It  is  interesting  to  look  at  the  solution  of  (3.3)  as  k 
increases  to  infinity,  thus  effectively  encountering  the  right  hand  branch.  We  anticipate 

l 

that  the  growth  rate  tends  to  a  constant  P£ ,  so  that  we  expand  the  growth  rate  as 


P  —  (po  +  k  *pi  -1 - ^ 


(3.4), 


the  size  of  the  second  order  term  is  found  from  the  size  of  the  layer  found  to  occur 
adajacent  to  t]  =  0,  from  elementary  observations  of  the  numerical  results. 

Thus  we  introduce  a  boundary  layer  variable  ,\'j  given  by 

7]  =  k~*x-  (3-5) 

In  this  layer  it  is  known  that  the  basic  flow  v  tends  to 


1  ■ —  A  k  3  x  4* 
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where  A  =  —  v9|q=o  =  Expand  the  disturbance  U  as 


U  —  Uq  +  k  3U\  -{-  •  •  •  . 


(3.7) 


at  leading  order  we  obtain 


~PoUq  +  2A  Uq  —  0. 


(3.8) 


This  is  a  compatibility  condition,  so  for  non-trivial  solutions  we  require  fio  =  2A,  at 
next  order  we  obtain 

W~(X(x  +  h))Uo  =  0\  (3-9) 

For  decaying  solutions  as  x  oo,  the  solution  of  (3.9)  is  given  by . 


U0  =  IC  Ai 


(3.10) 


To  solve  the  boundary  condition  at  x  =  0  we  require  that 


Ai 


(3.11) 


thus  we  can  find  the  second  order  growth  rate  Pi  in  terms  of  the  eigenvalues  of  the 
Airy  function.  Notice  that  this  function  has  many  zeros  on  the  negative  real  axis,  we 
are  interested  in  the  highest  value  of  the  growth  rate,  that  is  the  first  eigenvalue.  Thus 
as  &  — »  oo,  we  have 

V2A-^aU-§+...  ,  (3.12) 

33 

which  we  compare  with  the  right  hand  branch  growth  rate  in  section  4,  this  result  is 
also  included  in  figure  1,  note  this  agreement  is  reasonable,  but  it  should  be  noted  that 
£-3  term  does  not  tend  to  zero  very  quickly.  Also  included  in  figure  1  is  the  righthand 
growth  rate  value,  and  it  can  be  seen  that  the  solid  curve  is  tending  to  this  value. 


§4  Right  hand  branch  calculation 

We  consider  the  right  hand  branch  where  k  is  large,  and  large  Gortlcr  number, 
G  as  introduced  in  section  3.  Neutral  vortices  in  this  regime  were  described  by  Hall 
(1982  a).  In  the  neighbourhood  of  the  neutral  position,  temporally  amplified  modes 
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with  growth  rate  0{k 2),  are  found  to  occur.  In  this  regime,  we  know  for  neutral  modes 
G  ~  fc4,  so  we  write 

=  (4.1) 

We  note  that  the  vortex  activity  is  likely  to  be  concentrated  where  the  Rayleigh  cri¬ 
terion  is  most  violated,  for  the  problem  considered  here  this  occurs  at  the  wall,  which 
is  the  same  as  the  Taylor  problem,  one  of  the  problems  considered  by  Hall  (1982  a). 
We  consider  a  disturbance  proportional  to 

B  =  exp(^  +  i  /‘aW)  +  ‘A  (#)  +  •••<(#),  (4.2) 

where  1/e  =  fc,  and  thus  G  may  be  expanded  as 

OO 

G  =  £  9„£-4+".  (4.3) 


Thus  the  governing  equations  become 


(4.4a) 


(4.46) 


Here  we  consider  non-neutral  solutions,  thus  /?o  is  non-zero  it  is  possible  to  derive 
a  compatibility  relation,  so  in  (4.4)  we  now  have 

{U,V)  =  (uo,k~2v0)E,  (4.5) 

We  introduce  a  layer  of  thickness  A*  2G~»,  that  is  k~ *.  So  substituting  (4.5)  into 
(4.4)  in  this  layer,  located  at  the  cylinder,  we  obtain 

(At2+/?)At2u0  =  2vo,  (4.6a) 

(V2+/?)Vo  =  AuoAt2  (4.66) 

where  we  have  assumed  for  r)  <C  1,  v  — *  1  —  At?  +  •  •  Thus  combining  (4.6  a)  and  (4.6 
b)  we  obtain, 

(Af2  +  /?)  -  2A  =  0.  (4.7) 
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Thus  as  A*  — >  0,  the  growth  rate  j3  tends  to  %/2A,  which  agrees  with  the  inviscid 
problem  as  k  — >  co.  Note  that  we  also  have  an  order  one  cut  off  point  where  the 
growth  rate  is  zero  given  by 

Af  =  (2A)<.  (4.8) 

Thus  it  is  now  possible  to  continue  in  the  same  manner  as  Hall  (1982  a)  and  derive 
the  form  of  the  neutral  modes,  this  analysis  is  not  included  as  the  conversion  to  this 
case  is  trivial  and  has  the  same  conclusions.  Note  that  in  Denier,  Hall  &  Seddougui 
(1991)  they  found  that  the  position  of  the  vortex  activity  varies  with  A*,  we  do  not 
have  that  here  because  our  growth  rate  is  temporal  and  associated  with  8 /dr,  whereas 
in  Denier,  Hall  &  Seddougui  (1991)  the  growth  rate  is  associated  with  ud/dx.  They 
also  found  as  A*  — >  0,  that  there  is  a  singularity  in  the  growth  rate  and  thus  there  is 
an  intermediate  G*  wavenumber  regime  required  .o  match  with  the  inviscid  problem 
but  that  does  not  occur  here,  as  we  have  direct  correlation  between  large-  k  inviscid 
growth  rate  and  small- &  right  hand  branch  growth  rate. 


§5  The  numerical  work 


Now  let  us  discuss  the  numerical  scheme  employed  to  integrate  the  inhomogeneous 
system  for  U  and  V.  The  scheme  is  based  on  a  finite  difference  scheme  in  the  77-direction 
together  with  an  Euler  step  in  r.  Thus  for  example  (2.10  b)  becomes 


%ti-2  Vt  +  Vti 

h 2 


Vi 


=  V; 


*+l 


.2,7-  ,  VC 


2  h 


-  klV~  + 


Ar’ 


(5.1) 


where  i  denotes  grid  position  in  the  77-direction,  whilst  A  denotes  the  function  evaluated 
at  r  +  Ar  and  r  respectively.  Furthermore  h  and  Ar  are  the  step  lengths  in  the  77 
and  r  directions  respectively.  If  the  disturbance  velocity  field  is  known  at  time  r,  then 
the  tridiagonal  system  associated  with  (2.1)  can  be  solved  to  give  V  at  time  r  +  Ar, 
a  similar  discretization  of  (2.10  a)  then  produces  a  pentadiagonal  system  for  U,  this 
system  can  be  solved  for  U  at  time  r  +  Ar;  this  procedure  can  then  be  repeated  to 
march  forward  in  time.  This  scheme  is  fully  implicit  so  is  stable  for  Ar  =  0(h)  rather 
than  Ar  =  0(h 2)  which  would  be  the  case  with  an  explicit  scheme. 


12 


After  some  grid  checks  we  found  that  300  steps  in  f?,  with  infinity  r\  —  20  was 
sufficient  to  present  the  given  results  to  sufficient  accuracy.  The  value  of  Ar  used 
to  calculate  the  results  given,  was  Ar  =  0.001.  The  initial  value  of  r  where  the 
calculations  were  begun  was  varied  until  it  was  sufficient  small  to  calculate  U  and  V 
with  sufficient  accuracy,  the  value  used  was  Ar  =  0.01. 

As  an  indication  of  vortex  growth,  the  energy  was  computed,  defined  by 

OO 

E  =  J  U2drj ,  (5.2) 

o 


as  mentioned  in  the  introduction  and  then  the  growth  rate  was  calculated,  defined 

by 

<5-3> 


For  all  our  calculations  /?  was  found  initially  to  be  negative  so  at  sufficient  value  of 
time  the  flow  is  stable.  At  a  later  time,  say  r0,  we  found  that  /?  becomes  positive  until 
ri,  where  /?  again  becomes  negative.  Thus  there  is  a  finite  time  interval  where  the 
basic  unsteady  boundary  layer  is  unstable.  At  r  —  tq  and  r  =  t\  ,  the  instantaneous 
wavenumber  kr%  can  be  plotted  against  the  instantaneous  Gortler  number  Gri  to 
generate  the  neutral  curve  for  spin-up,  this  is  shown  in  figure  2,  for  G  —  1. 

We  conclude  that  spin-up  is  instantaneously  unstable  for  small  time  interval, 
~  i?e  3 ;  during  this  time  turbulent  flow  will  almost  certainly  be  set  up  so  that  stable 
rigid  body  rotation  may  never  be  achieved  in  the  practical  situation.  Furthermore 
we  note  that  the  induced  velocity  for  any  shaped  bump  can  be  found  by  numerical 
integration  of  the  transformed  velocity  field.  Several  different  humps  were  tried,  and 
the  results  were  found  to  be  quantitatively  similar,  so  we  only  include  the  plots  for 
the  humps  included  in  Denier,  Hall  &  Seddougui  (1991).  In  the  graphs  the  ordinate  is 
the  2-axis  and  the  abscissa  represents  the  time  axis.  Figure  3,  represents  the  contours 
of  maximum  V  component  for  the  hump  7r^ /4e— “2/16,  and  figure  4  the  contours. for 
Tr?z/32e~z2/16.  Note  these  are  basically  the  symmetric  and  anti -symmetric  cases. 
The  region  of  vortex  activity  is  basically  the  same  for  any  form  of  hump,  even  step 


functions. 
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§6  Conclusion 


We  have  shown  for  the  spin-up  problem,  that  Gortler  vortices  may  occur  when 
forced  by  an  axi-symmetric  hump  at  a  cylinder  started  to  move  impulsively,  also  that 
the  shape  of  the  hump  does  not  effect  the  amplitude  of  the  vortices  found  to  evolve, 
or  greatly  alter  the  time  of  the  onset  of  the  instability. 

We  note  that  the  growth  rate,  for  an  order  one  initial  wavelength  has  two  zeros, 
namely  t  =  tq  and  r  =  ti  ,  thus  between  To  and  t\  ,  (3  the  growth  rate  is  positive  and 
during  this  finite  time  interval  spin  up  flow  is  unstable.  The  neutral  curve  is  plotted 
for  the  humps,  notice  it  is  independent  of  hump  shape  as  it  is  done  before  the  Fourier 
inversion.  We  also  produce  figures  showing  vortex  growth  for  each  hump  considered, 
the  plots  are  of  time  against  the  maximum  value  of  V  for  a  given  z,  this  quantity  is 
taken  to  be  representative  of  the  vortex  amplitude.  Note  that  to  graphical  accuracy 
the  wavelength  of  the  vortices  obtained  corresponds  to  the  critical  wavelength  from 
the  relevant  neutral  curve. 

In  section  3,  we  described  the  inviscid  Gortler  modes,  which  we  found  to  have 
a  growth  rate  of  a/2A (7z,  as  k  — >  oo,  in  section  4  the  growth  rate  of  the  right  hand 
branch  in  the  limit  k  —y  0,  was  found  to  tend  to  G?V 2A,  thus  we  have  matching  of 
the  growth  rates,  there  is  therefore  no  further  structure  brought  on  by  a  singularity  in 
the  growth  rate,  as  was  found  to  occur  in  Denier,  Hall,  &  Seddougui  (1991). 

In  section  4  we  did  not  include  the  description  of  neutral  vortices  but  briefly  we 
see  that  the  first  order  disturbance  term  satisfies, 


d>v„  e* 

_  _  _y„  -  aVo  =  0. 


where 


and 


•-'•<'-*(1*1)) 


a  = 


3  a2 


where 


a  =  (  — 4  1  -  02  '  K°3 


'Y&  +  Y& 

Vot  V?u 


(6.1) 


(6.2a) 


(6.26,  c) 
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where  the  superscript  +  denotes  evaluation  at  rj+ ,  the  position  of  vortex  activity,  in 
this  case  equals  zero.  Thus  Vo  is  given  by  a  parabolic  cylinder  function,  ensuring 
the  required  decay  as  x  — ^ >  oo,  if  a  =  —  m,  thus  we  have 


r  2  da 

90  ~  Ql  —  92  —  ) 

Koo  voi  z  ^oo  Koi 

note  the  terms  V0'|"  are  the  coefficients  of  the  Taylor  series  of  the  basic  flow  evaluated 
in  a  layer  of  thickness  k~% ,  situated  at  r)+. 

It  is  also  possible  to  extend  the  small  wavelength  analysis  to  a  fully  non-linear 
calculation,  performing  a  similar  calculation  to  Hall  &  Lakin  (1988). 

It  should  be  also  noted  that  we  may  extend  this  analysis  to  further  geometries, 
the  first  obvious  extension  would  be  to  the  case  involving  a  sphere,  a  brief  outline  of 
the  method  involved  follows.  Firstly  it  is  necessary  to  invoke  slightly  different  scaling 
for  a  sphere  involved  in  spin-up  about  an  axis  of  symmetry,  the  scaling  are  now 

w  ~  W,  ~  aW2/v,  r  ~  a,  t~v/a2  p~pW2.  (6.4) 


These  scaling  are  the  same  as  those  employed  in  the  torsional  case  of  Otto  (1992  a), 
whereto  the  frequency  is  now  replaced  with  u/a2  Thus  the  governing  equations  become 


v  du 
r  d9 


r 


d 


r2  sin  9  d9 


(vsin  9)  ,  (6.5a) 


dv 


dv  v  dv  uv 


w 


1  dp 


2  du 


^-+R2Au^-  +  ~-—\-  —  cotO  =  --^+\V2v+~- 


dt 


dr  r  d9 


rdO 


r 2  d9  r2  sin  9 


(  dw  v  dw  uw  vw 


^  +  K‘Au^  +  ~  +  —  +  —  cotlH  =  |  V‘w- 


dt 


{ 


w 


dr  r  d9 


r 

1 


r2  sin2  9  ]  ’ 


d 


r  sin  9  d9 


(usin<?)  =  0, 


(6.56) 

(6.5c) 

(6.5d) 


where, 


V2 


i_d_( ad\  i  9  (  d\ 

r2  dr  y  dr  J  r2  sin  9  d9  y  d9  J 
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Note  that  we  have  assumed  that  the  flow  is  independent  of  the  azimuthal  coordinate 
< i i,  as  was  done  in  the  torsional  case  of  Ctto  (1992  a),  this  also  corresponds  to  assuming 
the  flow  in  the  cylinderical  case  is  independent  of  9.  A  Reynolds  number  has  also  been 
introduced,  where  this  is  given  by  Wajv. 

Prom  similar  arguments  to  those  employed  in  the  torsional  case  discussed  in  Otto 
(1992  a),  the  onset  of  a  vortex  stability  is  most  likely  to  occur  at  the  equator,  thus  we 
regard  the  basic  velocity  in  this  area.  The  derivation  of  the  basic  flow  was  attempted 
in  Benton  (1965),  and  he  showed  the  basic  flow  components  to  be  given  by, 


w  ~  sin  9f(r,t), 


(6.6  a) 


v  ~  cos  9  sin  9g(r,t). 


(6.6  b) 


Note  in  the  region  of  7r/2  we  introduce  a  boundary  layer  variable  9  <£  1,  and  the  basic 
flow  quantities  expand  as 

w  ~  f(r,t)  +  0(92) . . .  ,  (6.7a) 


The  Gortler  number  is  given  by 


ii  ~  0(9 )  +  •  •  • 


G  =  R2CTK 


(6.76) 


and  so  in  a  similar  manner  to  the  cylinder  problem,  .ve  employ  a  time  scaling  so  that 
the  Gortler  number  remains  of  order  one,  thus  we  introduce 

r  =  lJ-J7c5??,  t  =  i?c3r,  9=^  +  Rea9.  (6.9) 

Now  in  the  region  of  the  equator  of  the  sphere  the  basic  flow  is  found  to  be  given 
by  an  error  function,  we  perturb  this  basic  state  by  the  quantity 


(r7*u,Rc  *~',Rc  5w)  , 


(6.10) 


As  thus  we  arrive  at  the  governing  equations 


du  _  „  dp  ( d2  d2  \  „ 


(6.11a) 
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dv__dp  (&_ 

dr~  d9+  W  +  dd 2 )  V  ’ 

du)  „  /  d2  d2  \  „ 

Tfr +,JW~  \W  +  w) w’ 


(6.116) 
(6.11c) 
(6. lid) 


Thus  note  if  a  Fourier  transform  is  now  applied  to  the  these  equations  we  arrive  at  the 
equations  (2.9),  with  v  replaced  by  w.  Therefore  we  have  shown  that  a  sphere  involved 
in  spin-up  is  locally  unstable  to  the  same  vortex  structure  as  that  found  at  a  cylinder 
in  the  same  situation.  Although  it  should  be  noted  that  as  in  the  torsional  case  of 
Otto  (1992  a),  difficulties  may  be  experienced  actually  at  the  equator,  so  extension  of 
the  problem  into  the  complex  0-plane  may  be  necessary. 
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Figure  2  The 


neutral  curve 


Figure  2:  Neutral  curve  of  local  Gortler  number  Gt  against  local  wavenumbe 
ax,  for  a  form  of  hump  scaled  on  boundary  layer  thickness. 
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